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ABSTRACT 



In LTE downlink systems with a high frequency reuse factor, inter-cell interference has been viewed 
as a predominant factor limiting system performance. Therefore, we develop an iterative interference 
suppression receiver to mitigate inter-cell interference for LTE downlink systems in this paper. The key 
challenge here is how to accurately estimate the channel state information and spatial interference 
covariance matrix by using limited available pilots to improve the interference suppression capability and 
the successful probability of signal detection. To overcome the challenge, an iterative receiver making full 
use of the time-frequency correlation of channel is proposed. Simulation results show that the proposed 
iterative receiver can obtain a lower packet-error ratio compared with the existing algorithms. 

© 2013 Elsevier Inc. All rights reserved. 



1. Introduction 

In conventional cellular systems, the users at different locations 
from the same cell or from the different cells may employ the 
same frequency, which will inevitably result in intra-cell or inter- 
cell interference, respectively. The former can be generally avoided 
or mitigated by carefully designing transmission signals or allo- 
cating resource while the latter is hard to deal with if there is 
no cooperation among multiple cells [1]. Consequently, the inter- 
cell interference has been viewed as a predominant limiting factor 
for the performance of cellular systems, especially for the cellular 
systems with a high frequency reuse factor, such as long term evolu- 
tion (LTE) systems with full frequency reuse [1]. In this case, some 
effective schemes have been proposed to mitigate interference, 
which mainly include resource allocation and preprocessing at the 
transmitter side to avoid interference [1] and the interference sup- 
pression at the receiver side [2-11] to mitigate interference. In this 
paper, we mainly focus on the interference suppression at the re- 
ceiver side for LTE downlink. 

For the interference suppression receiver design, a basic idea 
is first to estimate the channel state information and interference 
covariance matrix based on available pilots, and then to imple- 
ment interference suppression and signal detection by using the 
estimated channel and the spatial covariance matrix [2-11]. This 
means that the interference suppression capability of receiver will 
strictly depend on the accuracy of channel state information and 
spatial interference covariance estimations. If the channel state in- 
formation is known, the receiver design only need to focusing on 
the signal detection, such as literatures [2-7]. Based on the known 
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channel information, the authors in [2] first use turbo principle to 
design a soft-input soft-output minimum mean-square error (MMSE) 
detector for code division multiple access (CDMA) systems while 
the authors in [3] use the iterative successive parallel arbitrated 
decision feedback (DF) to develop an MMSE-DF receiver for direct 
sequence CDMA systems. Furthermore, the authors in [4-6] ap- 
ply the turbo principle to MIMO systems and developed several 
efficient iterative MIMO receivers. Considering the time variation 
of the practical MIMO channel in a packet, a joint tracking and 
equalization algorithm based on a Kalman filter is proposed in [7]. 
Depending on the iterative process, the performance of the afore- 
mentioned receivers is continuously improved. 

However, the channel state information in practical systems is 
unknown and need to be estimated. Generally, it is difficult to 
obtain an accurate estimation of channel and spatial covariance 
matrix because the available pilots are limited and the interference 
is changed with the scheduled period in practical systems, such as 
in LTE systems [1[. To overcome this difficulty, several channel and 
spatial covariance estimation schemes which can be applied to LTE 
systems have been proposed in [8-11]. In [8], a time-domain low- 
pass weighting filter has been proposed to improve spatial corre- 
lation matrix estimation of interference, which is mainly based on 
the Cholesky decomposition to exploit channel frequency-domain 
correlation to smooth the noise. Unfortunately, this leads to the 
channel spectral leakage in practical systems so that the estima- 
tion performance is degraded. Therefore, a time-domain sine type 
weighing scheme with the frequency-domain moving average fil- 
ter is used in [9-11] to reduce the estimation error caused by 
the spectral leakage. The channel in practical LTE systems is gen- 
erally constituted by the time- and frequency-domain correlations, 
but the aforementioned schemes only use one part of both channel 
correlations. Therefore, we will propose an iterative interference 
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Interference cell 



Fig. 1. Multi-cell interference model. 

suppression receiver with enhanced channel estimation and im- 
proved spatial covariance matrix estimation in this paper. Its key 
idea is to make full use of the time-frequency correlations of 
channel to improve the estimation accuracy of channel state infor- 
mation and spatial covariance matrix in each iteration. Meanwhile, 
it can also use the iterative processing to continuously improve the 
performance of the proposed receiver. 

The rest of this paper is organized as follows. Section 2 intro- 
duces the signal model and the channel model in LTE downlink 
systems. In Section 3, an iterative receiver with enhanced channel 
estimation and spatial covariance matrix estimation is proposed. 
Simulation results are presented in Section 4 and conclusions are 
drawn in Section 5. 

Here are some notations used in this paper. 

• (-) H : the Hermitian operation; 

• E{-}: the expectation operation; 

• trace{-}: the trace operator; 

• diag{-}: the diagonal operation; 

• I m : the m x m identity matrix; 

• F m : the m x m Fourier transform matrix; 

• [A] n , m : the element at the nth row and the mth column of 
matrix A; 

• det(A): the determinant of matrix A; 

• 1 : /: from 1 to /. 

2. System model 

In this section, we first describe the signal model and then dis- 
cuss channel model and pilot structure in LTE downlink. 



the fcth subcarrier of the nth OFDM symbol, z„ j { is an N r x 1 noise 
and interference vector. The noise and interference vector, z nk , is 
assumed with zero mean and covariance matrix, R„ ,/<, at the fcth 
subcarrier of the nth OFDM symbol. The spatial covariance matrix 
of interference and noise, R„ k can be expressed as 



r CID 
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(2) 



where r, 



n,k 



denotes the correlation coefficient of the noise and 



interference from different received antennas, r\ and V2- It should 
be noted that we assume the transmit signal, the interference, and 
the noise are independent each other in this paper. 

Based on the Gaussian approximation [12], the probability den- 
sity function (PDF) of the received signal in (1) can be expressed as 
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From the above equation, the maximum likelihood detection (MLD) 
can be used to estimate the transmitted signal, x„ k, if the spatial 
covariance matrix, R„ ^, and the channel matrix, H„ k , are known. 
However, R„ ik and H,,^ in (3) are generally unknown and need to 
be estimated by using pilot signals or training symbols in practical 
systems. Therefore, joint estimation for x n k , R nk , and H„ /< is re- 
quired to maximize the PDF, p(y n ,fc|x n j ( , H„,/ ( , Rn,k). Since this joint 
estimation has very high complexity, we will develop an iterative 
scheme with low complexity in Section 3 to simplify it. 

2.2. Channel model and pilot structure 

To obtain a good result of the channel estimation and the spa- 
tial covariance matrix estimation, we will first discuss the channel 
model and the available pilot structure for LTE downlink in this 
section. 

The impulse response of equivalent baseband channel model 
from the n t th transmit antenna to the n,-th receive antenna can 
be described by 



h 



(n r n t ) 



,(n r n t ) 
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(4) 



2.3. Signal model and problem description 

As shown in Fig. 1, we consider a downlink multi-cell system 
consisting of N ce u s cells, which share the same available spectrum, 
i.e., the frequency reuse factor is one. Each cell has a base station 
(BS) with JV r receive antennas and U users each with N t transmit 
antennas. At the BS, the scheduler will select users and allocate 
them to a subset of resource blocks (RB's). For orthogonal frequency 
division multiple access (OFDMA) exploited in LTE donwlink, each 
user within the same cell will be orthogonally allocated a differ- 
ent subset of RB's. Therefore, the intra-cell interference is avoided, 
but the inter-cell interference still exists if there is no coordination 
among different cells. In this case, the received signal at the fcth 
subcarrier of the nth OFDM symbol can be expressed as 



Yn.fc = Hn,kXn,k + Zn.k, 



(1) 



where H„,k is an N r x N t channel matrix at the fcth subcarrier of 
the nth OFDM symbol, x„ /< is an N t x 1 transmit signal vector at 



where r, and a^'l"" are the delay and complex gain of the /th 
distinguishable path, respectively. Correspondingly, the frequency- 
domain response of the baseband channel in (4) can be ex- 
pressed as 



„(n r n t ) 
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(5) 



Since each distinguishable path consists of many unresolvable 
sub-paths, it is usually modeled as a wide-sense stationary (WSS) 
narrowband complex Gaussian processes with the same normal- 
ized correlation function. Different transmit-receive antenna pairs 
are assumed to have the same power-delay profile due to the 
small antenna space in centralized multi-antenna systems [14]. 
As a result, the time-frequency correlation of the channel can be 
expressed as [13] 



r H (Af, At) = a%r t (At)r f (Af), 



(6) 
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Fig. 2. Cell-specific pilot structure of two antenna port with normal cyclic prefix [15]. 
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Fig. 3. Iterative receiver with channel and spatial covariance estimation. 



where r t (At) is the time-domain correlation function determined 
by the Doppler spectrum, r/(A/) is the frequency-domain corre- 
lation function determined by power delay profile and is defined 

cr 2 , 

as r/(A/) = 2^1 T" ex P ( — j27r/T)), a, is the power of the /th dis- 

tinguishable path, er^ = J2i=o a f ' s tne t0,:a ' avera S e power of the 
time-domain channel, L is the length of the time-domain channel, 
i.e., the number of the distinguishable channel paths. 

Assume that there are /( subcarriers in an OFDM symbol and 
the subcarrier space is A/, then the duration of each OFDM 
symbol and the sampling interval are T s = 1/A/ and YKJ' re ~ 
spectively. Consequently, the equivalent discrete complex baseband 
channel can be expressed as 




where H n [ is an element from the n r th row and the n t th col- 
umn of Hn /; in (1). 

The distinguishable paths in practical systems are generally 
sparse [19,20], i.e., the time-domain channel has only limited non- 
zero paths. Therefore, we can exploit this fact to improve the 
estimation accuracy of the channel and spatial covariance matri- 
ces. 

Depending on the channel statistics in LTE systems, four dif- 
ferent pilot structures are defined for the downlink [15]. In this 
paper, we only exploit the cell-specific pilot structure, as shown 
in Fig. 2, to estimate the channel matrix and the spatial covari- 
ance matrix. For the cell-specific pilots, each RB only includes four 
pilots corresponding to each antenna port. Therefore, a key prob- 
lem is how to improve the channel and spatial covariance matrix 
estimation by using the limited pilots in each scheduled duration 
since different scheduled durations maybe face different inter-cell 
interference. 



3. Iterative receiver with interference suppression 

To maximize the PDF, p(y n ,k\Xn,k> H„,k» Rn,k). in (3), we need to 
jointly estimate the channel matrix, the spatial covariance matrix, 
and the transmit signal. However, the complexity of joint estima- 
tion is too high, we therefore develop a low-complexity iterative 
estimation scheme as shown in Fig. 3, which mainly includes three 
parts: the channel estimation, the spatial covariance matrix esti- 
mation, and the signal detection. In this scheme, we mainly use 
the time-frequency correlation of the channel to improve the accu- 
racy of the channel and spatial covariance matrix estimation, and 
then the enhanced channel estimation and the improved spatial 
covariance matrix estimation is used to improve the interference 
suppression capability and the successful probability of signal de- 
tection. 

3.1 Channel estimation 

In this section, we will develop a practical enhanced channel 
estimation scheme, which mainly uses the time-frequency corre- 
lation to improve the estimation accuracy. 

To the best of our knowledge, if the statistical information of 
the channel, the interference, and the noise is known, then MMSE 
estimator can be used to achieve the optimal estimation perfor- 
mance [13,14], However, the other statistical information except 
the noise is hardly obtained in practical systems. Therefore, least 
square (LS) estimator is generally used to obtain an initial esti- 
mation result, and then this result can be used to estimate the 
corresponding statistical parameters to further improve channel es- 
timation accuracy. 

Assume that K, p , K-d, and /C u are the index sets of the 
frequency-domain pilot subcarriers, data subcarriers, and unavail- 
able high frequency and direct current subcarriers for LTE down- 
link, respectively. Denote Af p and A/j to be the index sets of the 
OFDM symbols with and without pilots for LTE downlink, respec- 



90 



J. Fan et al. / Digital Signal Processing 24 (2014) 87-94 



tively. Then, LS estimation of the channel vector corresponding to 
the n t th transmit antenna can be expressed as [13,14] 
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(8) 



where k e tC p and n e A/p for the initial channel estimation stage 
while k e AC P U /Q and n 6 A/p U A/j for the iterative channel es- 
timation stage, H^" f * denotes the channel vector from the n t th 
transmit antenna to all receive antennas in the desired cell and 
can be viewed as the n t collum of H„ /<. In the following, the ini- 
tial channel estimation and the iterative channel estimation will 
be identically performed. For simplifying expressions, we denote 
the index set of all subcarriers corresponding to known available 
signals and the index set of all OFDM symbols corresponding to 
known available signals are defined as 



y~ A 

K, a — 



and 



M a 



AC P , for initial estimation, 

/CnU/Qj, for iterative estimation, 



A/p , for initial estimation, 

Nn U Mi, for iterative estimation, 



respectively, where the available known signals include pilots 
or/and estimation signals. 

After the channel state information at the known subcarri- 
ers, which are either the pilot subcarriers for the initial estima- 
tion or the pilot and data subcarrier for the iterative estimation, 
is obtained, the channel state information at the other subcarriers 
can be calculated by using the interpolation. The two-dimension 
MMSE channel interpolation for grid pilots has a very good perfor- 
mance [16], but it requires channel statistics and is very complex. 
Therefore, we will propose a robust and low-complexity chan- 
nel interpolation scheme consisting of the cascade of two one- 
dimension interpolations. 

Assume B is a frequency-domain interpolation matrix, then the 
channel vector on all subcarriers of the nth OFDM symbol from 
the n t th transmit antenna to the n r th receive antenna, H^J^ = 



[H ( n y d H ( n " r 2 n,) ■ ■ ■ H^">] , can be estimated by 
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(9) 



where AC = ]C P U ACd U ACu denotes the set of all subcarriers, while 
jj("r"t) denotes a channel column vector from all known subcarri- 

n,K, a 

ers, AC a . 

From (8) and (9), we can express the estimation mean-square 
error of the frequency-domain channel vector from the n t th trans- 
mit antenna to the n r th receive antenna at the nth OFDM sym- 
bol as 

= E{trace{«^ -B<^)«^ 
= E{trace{((Fx;,£ - BF KaX )h n nrn ' ) - Bz nXa ) 
x ((F K ,£ - BFK oX )h^ f) - Bz^y } } 



= trace} (F KX - BF KaX )C h (Fic X - BF KaX r 
+ BT nXa B H }, 



(10) 



where C h =E{h { " rn ' ) h n n ''"' ) "} contains the statistics of the time- 
domain channel from the n t th transmit antenna to the n r th receive 
antenna at the nth OFDM symbol and T n ,K. a = E{z„x 0 ZnX 0 H } ' s 
the frequency-domain statistic of interference and noise at all 



the known subcarriers. As mentioned before, different transmit- 
receive antenna pairs have the same channel statistics, then, 
the interference channels to the different receiver antennas will 
have the same statistics. We therefore have the same C;, for all an- 
tenna pairs and the same T n ,jc p at the different receive antennas. 

By minimizing the mean-square error in (10), we have the op- 
timal frequency-domain interpolation matrix as 



B = F KX C n F» aX {? KaX C h F» aX +T nXo ) 1 



^^(C-VFj^T^F^) 



r-1 



/Cq,/Z n,/Ca' 



(HI 



where c denotes a submatrix of Fourier matrix, F, with sub- 
scripts AC and C, C denotes a set of non-zero channel paths. 
It should be noted that the second equality uses the matrix in- 
verse theorem and we assume the correlation matrix of the time- 
domain channel, Q, and the correlation matrix of the interference 
and noise, T n< jc a are reversible. Theoretically, the interpolation ma- 
trix in (11) is optimal, but it is very hard to apply this scheme to 
practical LTE systems since obtaining the second-order statistical 
characteristic of the channel is difficult and a lot of high frequency 
subcarriers are unused in LTE systems. In fact, the frequency- 
domain correlation of the channel is mainly determined by the 
limited length of the time-domain channel [14]. Therefore, we will 
simplify the interpolation matrix in (11) as 



B = Fy xl 



K,C l C,K r K.,K 



u. 



(12) 



where U is a linear interpolation matrix, I . is a selection matrix 

with subscripts C and AC, and C is the index set of the selected 
non-zero channel coefficients. In fact, the interpolation matrix in 
(11) is mainly implemented by the following four steps: (1) the 
estimated channel at the known subcarriers is interpolated by us- 
ing the linear filter, U; (2) the frequency-domain channel based 
on linear interpolation is transformed into the time-domain one; 
(3) the non-zero time-domain channel coefficients are selected ac- 
cording to the set £; (4) the selected time-domain channel is 
transformed into the frequency-domain one. Obviously, the final 
interpolation result mainly depends on the selection of the non- 
zero time-domain channel coefficients, i.e., the determination of 
the set C. Since the estimation channel is a superposition of the 
perfect channel and the noise, the selected scheme can be deter- 
mined by 



£={/||F«,UH^ ) | 2 >P^ 2 , Kia) 



(13) 



where p is a proportional factor to select the non-zero channel 
paths and cr n 2 is the noise covariance. Essentially, this frequency- 
domain channel estimation scheme mainly uses the limited length 
of the time-domain channel to improve the channel estimation ac- 
curacy. 

In the above, we have obtained the frequency-domain channel 
estimation. In the following, we will develop a time-domain chan- 
nel interpolation scheme. Assume that the time-domain channel 
interpolation matrix is A, then the channel vector of all OFDM 
symbols at the fcth subcarrier from the n t th transmit antenna to 



the n r th receive antenna, H^- r ^ : 
be estimated by 
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where M denotes a set of all OFDM symbols in a transmission time 
interval and \M\ denotes its cardinality, H?Jf n ^ denotes a column 
channel vector from all known symbols. Similarly, by using the cri- 
terion of the mean square-error minimization, we can obtain the 
time-domain interpolation matrix as 
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(15) 



where C 
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ii7/ H (nrnt) H (n r nt) f 



corresponds the channel time- 



varying characteristic, Cj = E{z 



i(nrn f ) ; (n r n t ) r 
■Af,k L M,k 



} reflects the statistic 



of frequency-domain channel estimation error, lj\f a ,j\f ' s a selection 
matrix with \M a ,M^HaM H = Based on (15), the channel at the 
subcarrier k corresponding to the different OFDM symbols will be 
obtained. 

As mentioned before, it is difficult to get second-order statistics 
of the channel, we will therefore propose a robust time-domain 
channel estimation scheme. Assume that the carrier frequency 
offset due to the oscillator mismatch between the transmit- 
ters and the receivers has been compensated, then the channel 
time-varying characteristics will be mainly determined by the 
Doppler spectrum. Consequently, we can approximately express 
C H as [13] 



c H [0] 
c H [l] 



Cff[l] 

c H [0] 



.c H [W\-1] c H [\Af\-2] 



c H [W\-l] 
c H [W\-2] 

c H [0] . 



(16) 



where c H [n] = sin n ( ™ d) , do d = f d is the maximal Doppler off- 
set determined by the mobile speed. Because the mobile speed 
of each user can be estimated easily, we assume that the user 
speed is known. Meanwhile, for simplicity, we can approximate 
C? as 



C 2 «a n 2 I. 



(17) 



Obviously, the mobile speed of each user and the noise covari- 
ance can be easily estimated, therefore, this scheme can be easily 
applied to practical LTE systems. 

3.2. Spatial covariance matrix estimation 

In the last section, the channel estimation accuracy is en- 
hanced by using the time-frequency correlation characteristics of 
the channel. Can we use this characteristics to improve the ac- 
curacy of the spatial covariance matrix estimation? To answer this 
question, we will first analyze the time-frequency correlation char- 
acteristics of the spatial covariance matrix. If the answer is yes, we 
can further use it to improve the accuracy of the spatial covariance 
matrix estimation. 

3.2.1. Analysis of the spatial covariance matrix 

Without loss of generality, we assume after CP of the received 
signal is removed, the frequency-domain expression of the inter- 
ference from the other cells and noise at the n t th receive antenna 
can be written as 



»("r) 
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„("r) 



(18) 



where v„ is a noise vector with zero-mean and covariance matrix 
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(19) 



is an asynchronous time-domain transmit signal from the n\ inter- 
ference cell, 
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(20) 



is the channel matrix from the BS of the n/th interference cell to 
the n r th receive antenna of the desired user, K, K g , m, and N; are 
the number of the subcarriers in an OFDM symbols, the CP length, 
the asynchronous time-offset, and the number of all interference 
cells, respectively. Obviously, if m = K, then the interference signal 
will synchronize with the desired signal. 

From (18), we can express the channel correlation of arbitrary 
two receive antennas as 



l n,K 
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In practical LTE systems, the transmit signals from the different 
cells are independent each other, we therefore have 



E{s<' 



("nn/i) c (nr2";2r 

3n 



I, n n =n/ 2 , 
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(22) 



Then from (21) and (22), we can simplify the channel correlation 
of arbitrary two antennas as 



Jnr\n r2 ) 
l n,K 
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n,=l 
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(23) 



From (23), we find that the element, r^J ( 3." r2 \ of the spatial covari- 
ance matrix is similar as the channel correlation. To show this fact 
more clearly, we give an example with K = 512, K g = 36, m = 256, 
INR = 10 dB in Fig. 4, where INR denotes the interference-to- 
noise ratio. From this figure, we find that time-domain length of 
spatial interference covariance matrix is limited which is similar 
as the characteristic of time-domain channel no matter whether 
the interference are synchronous with the desired signals or not. 
It should be noted that in this figure, the frequency-domain zero- 
zone is because the high frequency subcarriers are unused while 
the time-domain zero-zone is because the time-domain channel 
has limited length. These analysis shows that we can also use the 
time-frequency correlation to further improve the accuracy of spa- 
tial covariance matrix estimation. 

3.2.2. Improved spatial covariance matrix estimation 

Once the channel information and the transmitted signal are 
given, the instantaneous spatial covariance matrix at the fcth sub- 
carrier of the nth OFDM symbol can be estimated by 



Rji,/< = (Wn,k - Hn,kXn,k)(yn,k ~ ^n.k^n.k) 



(24) 



92 



J. Fan et at / Digital Signal Processing 24 (2014) 87-94 



0.5 - 




8 0.1 



£ 0.5 



O 



- Synchronization 
Asynch ronization 




200 300 
Sub-carrier Index 



400 



500 



- Synchronization 
Asynchronization 



100 



200 300 
Time Index 



400 



500 



Fig. 4. The time-frequency correlation of spatial interference covariance matrix for 
LTE downlink with K = 512, K. = 36, m = 256, 1NR = 10 dB. 



»(n ri nr 2 ) _ „(n r2 n ri )H 
T n,K ~ T n,K 



(29) 



where f, 



("r 2 n ri ) 
n.K 



can be obtained according to (28). Then from (28) 



and (29), the corresponding spatial covariance matrix can be esti- 
mated by 



Vi,/c : 



r n,k 



? (lN r ) -, 
'n,k 



? (N r JV r ) 



+ ah. 



(30) 



'n,k 'n,k - 1 

It should be noted that the filter matrices, A and B, here have 
the same design methods as the channel estimation, but are with 
different filter parameters. 

3.3. Signal detection 

Once the estimation results of channel state information and 
spatial covariance matrix are obtained, we can estimate the trans- 
mitted signal by maximizing the condition PDF in (3). Conse- 
quently, the estimated signal can be expressed as 



Based on the instantaneous estimation in (24), several schemes 
in [8-11] have been proposed to obtain the spatial covariance ma- 
trix estimation. 

In [8], the authors assume that all OFDM symbols have the 
same instantaneous spatial covariance matrix so that the statistical 
spatial covariance matrix can be estimated by 



1 N 
JV 



k- 



(25) 



Based on the estimation results in (25), the Cholesky decomposi- 
tion can be used to further improve the accuracy of the spatial 
covariance matrix of interference [8], If the instantaneous spatial 
covariance matrices are different at different OFDM symbols but 
are very close at several adjacent subcarriers [9-11], the statistical 
spatial covariance matrix can be estimated by 



*n,k : 



Q/2 
a=-Q/2 



(26) 



where a q = i [9-11]. Then, the estimation results in (26) can 
be further smoothed to improve the estimation accuracy by using 
the frequency-domain correlation of channel. However, the afore- 
mentioned schemes only use one part of the time-domain correla- 
tion of channel and the frequency-domain correlation of channel. 
Therefore, we will simultaneously use the time-frequency correla- 
tion of channel in the iterative processing to improve the perfor- 
mance of receiver. 

As mentioned before, the correlation characteristic of each 
transmit-receive antenna pair has a limited length. Then, we can 
use the similar idea as the channel estimation to improve the spa- 
tial covariance matrix estimation. Correspondingly, we can express 
the frequency-domain estimation as 



iv, ) 



-(n ri nr 2 )_ j-OVrj 
l n,K — Dr n.K a 



(27) 



where f, 



n,K„ 



is the correlation characteristic of the n r . th receive 



antenna and the n r2 th receive antenna at the known subcarriers 



K, a . Meanwhile, the time-domain estimation can be written by 



~(n rl n r2 ) _ fl ~(n r in r2 ) 
r M,k ~ M a ,k > 

for n r i ^ n,-2. While for n rl < n r 2, we have 



(28) 



x n , k = min |(R n !. /2 ) H (yn,i<-H nik x n>/< )| 



(31) 



where X denotes a set consisting of all possible transmitted 

-1/2 H 

signals, x n ,k- In fact, the operation of (R n / ) in (31) is to whiten 
the interference while improving the successful probability of sig- 
nal detection. Therefore, we need to continuously improve the ac- 
curacy of the spatial covariance matrix in the iterative processes. 
Correspondingly, we can also express the bit log-likelihood ratio 
(LLR) of the estimated signal as 



Urn) 
L n,k 



ln- 



E. 



^(x„, fc |y n ,/ ( .R n j/ 2 ,H nT , ( ) 



E. 



Xn.ke^m,- 



(32) 



where X m ±\ denote the sets of all possible bit vectors with the 
mth element as ±1 and 



When performing iteration, we can use the initial LLR estima- 
tion, L^™\ as an a priori information to update a posteriori LLR 



information, L 



D,(m) 



fD,(m) _ f 
L n,k - L 



n.k 

A x (m) 
n.k 



which can be expressed as [17,18] 



+ ln 



Ex n . 
Ex„ 



e ^(x„,i i ,x„_( [ _| m ],L°j [m] |y„ i i.,R I1 j /2 .H„_, ( ) 
s HXn,k,Xn,k,lm]X°x lm] \yn,k,\l /2 ,Hn,k) 



(33) 



;El(m) 



where 



^(Xn,k> Xn,k,[m], L n k [ m ] j y n ,k, , Hn.k) 

= -\( R n,k ) (yn,k-H n , k Xn,k)| + -x' n 



/<,[m] L n,/<,[m]' 



^n /< m ' ' s tne ext:rms i c information, x„ ;< r m ] is the subvector of 
x n,k which is obtained by removing the mth element of x n k, 
and L n 1 fc j . has the same form as X n ,k,[m]- By the iterative pro- 
cessing, a posteriori probability of the received signals can be con- 
tinuously maximized. 
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Table 1 

Simulation parameters. 



Carrier frequency: 2 GHz 


Sampling frequency: 7.68 MHz 


TO duration: 1 ms 


OFDM symbols per frame: 14 


FFT size: 512 


Subcarrier spacing: 15 kHz 


System bandwidth: 5 MHz 


RB: 12 subcarriers x 7 OFDM symbols 


Code: Turbo 


Received antenna: 2 


Modulation 


Desired cell: 16QAM 




Interference cell: 64QAM 


Transmit data streams 


Desired: 1 




Interference: one is strong and others are weak 


Channel model (19,20] 


Type: ETU 




Speed: 30 km/h 


Proportional factor: p 
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Fig. 5. The packet error ratio versus SNR when SIR = 4 dB. 

4. Simulation results 
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Fig. 6. The packet error ratio versus SIR when SNR = 5 dB. 
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In this section, we compare the performance of our scheme 
with that in [8-11], For simplicity, we call the scheme in [8] and 
the scheme in [9-11] as Li's scheme and Jang's scheme, respec- 
tively. In our simulations, we exploit the three-cell system models 
and assume there is a strong one in the interferences from the 
adjacent cells while others are weak. Denote P s and Pj as the 
power of transmit signal and the power of the strong interference, 
respectively, then we can define signal-to-interference ratio (SIR) as 
S1R= ^. Meanwhile, 5 MHz system bandwidth, Turbo code, and 
QAM modulation are exploited. The detailed simulation parame- 
ters are listed in Table 1. 

In Figs. 5 and 6, we show the packet-eiror rate (PER) versus the 
signal-to-noise ratio (SNR) when SIR is 4 dB and the PER versus 
SIR when SNR is 5 dB, respectively. From these figures, the pro- 
posed scheme has a better PER performance compared with Li's 
scheme and Jang's scheme. Meanwhile, the initial spatial covari- 
ance matrix in Li's scheme is estimated by the proposed scheme, 
this therefore makes Li's scheme to have a better performance 
than that of Jang's scheme. This also further verifies that the pro- 
posed scheme outperforms Jang's scheme. It should be noted that 
although the proposed algorithm may use the iterative process 
to continuously improve the estimation accuracy of the channel 
and the spatial covariance matrix, it will pay a price of increas- 
ing computation complexity compared with Li's scheme and Jang's 
scheme. 

In Fig. 7, we show the PER performance versus the number of 
iterations when SIR is 4 dB. From this figure, the PER performance 
of the proposed scheme will decrease as the number of itera- 



Fig. 7. The packet error ratio versus the number of iterations when SNR = 5 dB. 

tions increases. Meanwhile, the PER performance of the proposed 
scheme with the high number of iterations will rapidly decrease as 
the SNR increases. 

5. Conclusions 

In this paper, we have proposed an iterative receiver with inter- 
ference suppression for LTE downlink systems. It makes full use of 
the time-frequency correlation of channel to improve the estima- 
tion accuracy of channel and spatial covariance matrix. With the 
more accurate channel and spatial covariance matrix estimation, 
the capability of interference suppression and the successful prob- 
ability of signal detection can be further improved by the maxi- 
mum a posteriori detection. Since the proposed scheme is based 
on the pilot structure for LTE downlink systems and considers the 
practical constraints, it can be easily applied to receivers for LTE 
systems. 
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